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Abstract 

This paper presents a new connection between the generalized Marcum-Q function and the confluent 
hypergeometric function of two variables, $3. This result is then applied to the closed-form charac- 
terization of the bivariate Nakagami-m distribution and of the distribution of the minimum eigenvalue 
of correlated non-central Wishart matrices, both important in communication theory. New expressions 
for the corresponding cumulative distributions are obtained and a number of communication-theoretic 
problems involving them are pointed out. 

Index Terms 

Marcum-Q function, confluent hypergeometric functions, bivariate Nakagami-m, non-central Wishart 
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I. Introduction 

The number of special functions that make appearances in the communication theory arena continues 
to grow. Some of these are tabulated and well studied, with readily available approximations, bounds, 
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and asymptotic expansions. Other special functions, however, are not so well studied and their numerical 
computation is not always as accurate and efficient as would be desirable. In such cases, the establishment 
of connections with other special functions can greatly reinforce the analytical toolbox available to 
researchers. 

A wealth of works have studied the Marcum-Q function (see, e.g., |T]]-||5j), deriving useful bounds 
and approximations and evidencing applications thereof. Routines for its accurate and efficient evaluation 
have been extensively developed and, in fact, the generalized Marcum-Q function is included in most 
common mathematical software packages. 

Confluent hypergeometric functions also appear in a fair number of problems within communication 
theory and signal processing (6}-(9j pOj . Moreover, recent distributional results J8j, fTT| show that a 



number of other such problems (see, e.g., j7), p2j-[22J) can be solved in terms of confluent hypergeo- 



metric functions. Chief among these stands the $3 confluent hypergeometric function |23|, |24|, which 
does not lend itself to easy and precise evaluation. 

This paper shows that <3?3 can be expressed in terms of the generalized Marcum-Q function. This con- 
nection is then applied to the closed- form characterization of two important distributions in communication- 
theory: the bivariate Nakagami-m distribution, and the distribution of the minimum eigenvalue of cor- 
related non-central Wishart matrices. New expressions for the corresponding cumulative distribution 
functions (CDFs) are obtained and some communication-theoretic problems involving them are pointed 
out. In particular, the connection unveiled herein, in combination with recent results in pT) , settles the 



standing conjecture p3| p4] p. 174] that the bivariate Nakagami-m CDF can be expressed in terms of 



the generalized Marcum-Q function. 



II. Preliminaries 



This section introduces the special functions under study and recalls a few properties of interest for 
the derivations that follow. 



Definition 1: The generalized Marcum-Q function is defined as |14|, [25 1, |26| 



Q m (a,b) = J^ -^ri ex P (~ a 2 X 1 Im-l {ax) dx (1) 

where a > and b > are real parameters and I m {-) is the mth order modified Bessel function of the 
first kind. The order index m is an integer and typically m > 0, yet ([T]) holds for negative orders too and 
a useful relationship between Marcum-Q functions with positive and negative orders has been reported 
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in (27}, namely 

Q m (a, b) = 1 - Qi- m (6, a) . (2) 



Definition 2: The $3 confluent hypergeometric function of two variables is denned as [23 Eq. 9.261.3] 

* s (6 iqw ,*) = ££^^. (3 ) 

where 6,c,to,z G I, 0, —1, —2, and (t) r = F p^ r ^ denotes the Pochhammer symbol with T(-) the 
Gamma function. The $3 function is one of the bivariate forms of the confluent hypergeometric function 



i-Fi(-, •;■) 1 23 1- Note that $3 does not exist for non-positive integer values of c due to the singularities 
of the Gamma function. Next, we introduce a regularized version of this function, which is valid for any 

cel. 

Definition 3: The regularized <3?3 function is defined as 

1 

$3 (b, c; w, z) = ^7^3 (b, c; w, z) 

r(c) 

= y* (% ^ m (4) 

2-^2-^> r (c + k + m) k\m\ 

fc=0m=0 V y 

with 6, c,w, z £ E. For the special cases it; = and/or z = 0, $3 reduces to 

l> 3 (6,c;0,z) = z( 1 - c )/ 2 / c _ 1 (2^) (5) 
$ 3 (6,c;u;,0) = iFi(6,c;u;) (6) 
$ 3 (6,c;0,0) = -^ (7) 

r c 



where c;w) = iF\(b, c;w)/T(c) is the regularized confluent hypergeometric function [23 1. 

The Laplace transforms of $3 and $3 are known and, due to their simple form, will be crucial in 
the ensuing derivations. Given the function f(t) = i c_1 <l>3 (6, c; xt, yt), its Laplace transform is given by 
(24) Eq. 4.24.9] 

C{f(t); S }=s- c (l-^y b > s . (8) 
III. Main Result 

The main result, presented in this section, rests on two new lemmas that provide, respectively, a new 
representation for the generalized Marcum-Q function and a recursive relationship for 63. 
Lemma 1: The generalized Marcum-Q function can be expressed in terms of $3 as 

Q m{ a,b)=( a ^ 1 ^(-^^Wl^-m;^,^), m G Z. (9) 
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Proof: See Appendix [A] 
Lemma 2: The regularized $3 function can be obtained recursively as 



ZN^fc? 1 



§z{b,c;w,z) = (-) V -A(fe,c;z)4» 3 (l,c-i;w,z) (10) 

j=0 



for any 6 G Z and b > 0, with «4j(6, c; z) being the polynomial on z given by 

A(h x (-l)^-i L ^ J (-l) fc (6 - j + fc) t _ fc ( C - i - 1 + fc),_ 2fc fc 

Proof: See Appendix [B] ■ 
Leveraging the foregoing lemmas, the main result in this paper can be put forth. 
Theorem 1: The $3 function is given in terms of the Marcum-Q function as 

i=0 v 7 



with 5, c G Z, /> > 0, z / 0, w / 0, and A.i(b, c; z) given by gl). 

Proof: The result follows directly from lemmas 111 and after solving for $3 (1, c — w, z) in d9 



and substituting in (10 1 



A few comments on the above theorem are in order: 

• For the special cases w = and z = 0, simpler connections are respectively given in ([5]) and ([6]) in 
terms of I m (-) and i-Fi(-, •; •). If u; = z = 0, then ([7]) gives $3 in terms of the Gamma function. 

• Albeit <l>3 (b, c;w, z) is defined for any real value of its arguments, Thm. |T| is restricted to integer 
values of b and c; these are precisely the cases of interest in communication theory. Moreover, the 



validity of ( [12] ) can be straightforwardly extended to c € R by applying ( |34] > in Appendix [A] in place 
of ©. 

• Negative values of to, 2 (which do not correspond to known communication theory problems) imply 
complex arguments of the Marcum-Q function. This is not an issue since the Marcum-Q definition 
in ([I]) also holds for complex arguments by analytic continuation |T|. 
Thm. [TJ allows expressing any result involving $3 in terms of the generalized Marcum-Q function. 
Besides having archival value, this relationship greatly facilitates both the evaluation of such results, and 
any subsequent analysis thereof. 

IV. Applications 

The generalized Marcum-Q and the $3 functions appear in a number of communication theory 
problems. The new connection between these functions presented in Thm. [T] can be therefore applied 
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directly to such problems. For instance, $3 appears in two important distributions: the bivariate Nakagami- 
m distribution, and the distribution of the minimum eigenvalue of non-central Wishart matrices. We next 
exemplify the applicability of Thm. [T] to these specific problems. 

A. Bivariate Nakagami-m Distribution 

Thm. [T] together with recent results for the bivariate Nakagami-m distribution in pl| , confirm the 



conjecture made in [15| and in [14 p. 174]: the joint CDF of two correlated Nakagami-m variables 
can be expressed in terms of the generalized Marcum-Q function, thereby completing the landscape of 
(closed-form) bivariate characterizations of the most common fading distributions. 

Corollary 1: Let Ri and R2 be two correlated Nakagami-m random variables with positive integer 
fading index m, respective variances Q\ = E[i2f] and ^2 = E[i?|], and correlation coefficient p, and let 
Ri and R2 be the normalized (unit-variance) versions of R\ and i?2- Then, the joint CDF of R\ and R2 
is given by 

P / \ l(m,mrl) - ror ; (mrf) fc 
V 7 ' fc=0 



m-lm-fe2(i-l) r - 2 \ fc + i -' r - 1 / 1 u\%\ 

fe=0 i=l r=0 \ r / \ I 



(13) 



where a = y jz^, b = a^fp, 7(-, •) is the lower incomplete gamma function |23|, and A r {-, •; •) are the 
polynomials defined in ( [TT] ). It follows that 

F KuKA r 1 ,r 2) = F^(^,^). (14) 

Proof: The bivariate Nakagami-m CDF is given in [11, Eq. 14] in terms of the $3 function. Then, 
( [T3] ) is obtained by virtue of Thm. [T] after normalization of the random variables and some algebraic 
manipulations. ■ 
Corollary 2: The well-known expression for the bivariate Rayleigh CDF can be recovered from Corol- 
lary [T] by setting m = 1 and using [23, Eq. 8.352.1] to expand the incomplete Gamma function, yielding 



F R u R 2 ( ri ' r ^ = 1 ~ e T2 ~ e ri Qi{ r 2a,nb) + e r2 Qi{r 2 b,r ia ) (15) 

consistently with [28, Appendix A]. 

The closed-form characterization of the bivariate Nakagami-m distribution had remained an open 



problem for decades. Existing expressions involved infinite summations [15]-|17| or an integral of the 
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product of Marcum-Q functions fT8~[ . The expression that recently appeared in pT[ , in terms of $3, has 
been rewritten by virtue of Thm. [T] into the convenient form given in Corollary [1] The problems whose 



analysis can benefit from this form include [11 1, 1 14|-fl22 



• Determining the impact of fading correlation in dual-diversity reception or transmission pT| , p9| 
sect. 5.2.5]. 



Analyzing the level crossing rate and average fade duration of sampled fading envelopes 1 11 1. 



Establishing the transition probabilities for a first-order Markov chain that models a fading process 



[21 1, [22]. This, in turn, can be applied to approximate the envelope of channels with non-independent 



fading [22| or to model the decoding success/failure with automatic repeat-request (ARQ) over 



successive channel realizations |21|. 



B. Minimum Eigenvalue Distribution of Non-Central Wishart Matrices 

Definition 4 (Non-central Wishart matrix): Let X be an n x m (n > m) random matrix distributed as 
CMn^mi"^ ') In ® E), where S is the covariance of the independent complex Gaussian row vectors of X, 
and T G c nxm _ Then W = X^X is a complex non-central Wishart matrix that follows the distribution 
CW m (n, S, 0) with = ^~ 1 X H T the non-centrality parameter. 

Concerning the extreme eigenvalues of W, distributional results are available for uncorrected central 
(S = I m , T = 0), correlated central (T = 0), and uncorrelated non-central (X = I m ) complex 



Wishart matrices (see, e.g., [29]-[33]). However, tractable results for the correlated non-central case had 
been unavailable until [8], where the minimum eigenvalue distribution of correlated non-central Wishart 
matrices has been expressed in terms of $3 for some special cases including a square X, i.e., for m = n. 
Specifically, the CDF of A m i n , the minimum eigenvalue of W ~ CW m (m, S, 5] _1 T^T), is given when 
T has rank one as (HJ 

F Amin (A) = 1 - exp (-7/ - Atr(S- 1 )) $3 {m,m, V , A/i) (16) 

where r\ = tr(0) and \i = tr(0S _1 ). 

Analogous expressions to ( [T6| ), also in terms of $3, are found for other special cases such as 2 x 2 
Wishart matrices with arbitrary degrees of freedom, i.e., m = 2 with arbitrary n, or 3 x 3 Wishart matrices 
with n = 4 degrees of freedom [8]. 

Using (16) and Thm. [T] we can express the minimum eigenvalue distribution of correlated non-central 



Wishart matrices as follows. 
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Corollary 3: The CDF of the minimum eigenvalue of W ~ CW m (m, S, S 1 y H T) with rank-one 
T can be expressed as 

F Amin (A) = 1 - exp (-r, - AtrOS" 1 )) T(m) (^j™ ' 

x A(m,m;A/i) exp ^ + — j Q 2 -m+» ( \/2r?, y 2 ~ j ( 17 ) 



where the polynomials Ar(-,-;-) are as in ( 11 ). 

Corollary [3] is restricted to the case of a rank-one non-centrality parameter, which however is typically 
assumed in multiple-input multiple-output (MIMO) communication systems with a direct line-of-sight 



path between the transmitter and the receiver |34|, |35|. Given the complexity of the underlying joint 



eigenvalue distribution, the CDF in ( 17 1 is remarkably simple, involving only a finite sum of generalized 
Marcum-Q and elementary functions. 

The minimum eigenvalue distribution is important in the analysis of MIMO channels J7J, p2j , [ 13 1 
where the received signal vector is modeled as 

y = Hx + n (18) 

with H G £_N r xN t j-jje channel matrix containing the gains between the N t transmit and 7Y r receive 
antennas, n G c^xi t ^ e no j se vector, and x G S the transmitted signal vector with entries drawn from 
an alphabet S. The minimum eigenvalue of H^H determines the minimum distance, d m - m , between the 
noiseless received signal vectors and, thereby, the error probability of a MIMO maximum likelihood 
(ML) receiver. It can be shown that p3| 



dmin > V Amin ^0 (19) 

where A m i n is the minimum eigenvalue of H^H and do is the minimum distance between the elements 
of S. Altogether, the performance of the MIMO ML receiver is strongly linked to the distribution of 
A m in> which is given in ( fT7| ) for H having non-central correlated Gaussian entries; this encompasses both 
Rayleigh and Rice fading with spatial correlation. Since ML becomes computationally unwieldy as the 
number of antennas or the transmission alphabet cardinality grows, linear and successive cancellation 



receivers become attractive. The performance of such receivers also depends on A m i n [12], |36|. In fact, 
the post-receiver SINR of the zero-forcing (ZF), minimum mean square error (MMSE), and Vertical Bell 
Labs Layered Space-Time (V-BLAST) receivers satisfies (12} 

SINR > Amin (20) 
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where E s is the energy per symbol, i.e., E[xx^] = (E s /N t )lN t , and a 2 is the noise variance. 

Combining Corollary [3] with (20), the outage probability of MIMO receivers can be analyzed in fairly 



broad generality. In addition, Corollary [3] has further applications, e.g., in the design and analysis of 
adaptive MIMO multiplexing-diversity switching [7] or, in the context of econometrics, in characterizing 
the weak instrument asymptotic distribution of the Cragg-Donald statistic (37| . 
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Appendix A 
Proof of LemmaQ] 



The generalized Marcum-Q function can be obtained as the contour integral in the complex plane [38 1 



Qm (a, b) = exp ("^-) f pm{ l_ p) Q (j + #p) ) dp (21) 

S J 

Z m (a,6) 

where <f Tg = f Tg and To is any closed contour enclosing the singularity at p = (in a counter- 
clockwise direction) and no other singularities of the integrand (cf. Fig. [TJ. For convenience, we express 
the integral X m (a,b) in (21 1 as 



l m (a, b) 



jf F(p) exp (^j dp (22) 



F(p) = -^—p- m e W (^). (23) 



with 

/•'•:/') , , 

1 — p \2p 

Let us consider the contour C depicted in Fig. [T] where c is chosen to the right of all the singularities 
of F(p). Letting R — > oo, the contour integral along C equals 

£ F(p) exp {^j dp jf +3 F{p) exp (^j dp + £ F(p) exp (^j dp. (24) 
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Fig. 1. Contour integration for integral X, 



Alternatively, we can apply the Cauchy-Goursat theorem to obtain 

b 2 p\ I w (b 2 p^ 



j> F(p)exp (jy^ d P = j> f (p) ex P (j~2~j dp + j> F(p) exp (^rj d P ( 25 ) 



where To and Y\ are closed contours enclosing the singularities at p = and p = 1, respectively. 

'by 



Combining d24b and (|2"5|), 



l 

27TJ 



r (p) exp 

C— JOO 



dp+ <j> F(p) exp ^^y^ dp = <j> F(p) exp ^^y^ ( 26 ) 



The first integral in (26) is related to the inverse Laplace transform of F(p) as 



1 

2wj 



F{p) exp 



C— JOO 



dp= £- L {F(p);t}\ 



The integral along Cp can be shown to be zero as follows. The modulus of F{p) in Cp is 

1 1 



|F(p)| 



p=ReJ e 



\1 — p\ \p\ 



J/P 



with 5 = a 2 / 2. Then, we can use the inequalities 

1 



1 — p 

S/p 



< 



1 



1 2 



,Re(<S/p) 



< e W*l < e 



|<5/-Ro| 



(27) 

(28) 

(29) 
(30) 
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with arbitrary (finite) Rq, to write 



K r-i 

Thus, \F(p)\ Reje < KR~ l for some I > on Cp and, according to Jordan's lemma, the integral § ( 



a. 



l m (a,b) = £- L {F(p);t} 



Res < F[p) exp 



(32) 



equals as R — > oo. Plugging (22) and (27 1 into (26 1 and applying the residue theorem, we arrive at 

where Res{-} p= £ denotes the residue at p = £. Then, we calculate the residue and use (|8j) to solve the 
inverse Laplace transform, which yields 



. , ,'& 2 \ m ~ / b 2 a 2 b 2 \ fa 2 + b 2 

Im(a,b) = - \ — \ $3 ( l,m + l;—,-j- J +expf — — 



(33) 



Further substituting (33 1 in (21 1 gives the alternative form for the generalized Marcum-Q function, 

,2\ m / a 2 + b 2 



Q m (a, b) = 1 



exp 



$ 3 l,m + 1; 



b 2 a 2 b 2 



, 2 J ""' \ 2 7 u V '"" ' ' 2 ' 4 

which is valid for any m 6 M. 

Finally, (|9]) is obtained by combining (34) and Q, completing the proof. 



(34) 



Appendix B 
Proof of Lemma[2] 

The recursive relationship in ( 10) is derived by using ([8]) and the frequency differentiation property of 
the Laplace transform, 

dF(s) 



£{t-f(t)} 



ds 



(35) 



Taking the first derivative of ([8]), in light of (35), 

i c l>3 (b,c;xt,yt) = 



C- 1 \ey' s [ hxs -(c+2) ( ± 



-(6+1) 



+ cs 



-(c+1) 



+ ys 



<c+2) ( 1 



(36) 



where C {•} stands for the inverse Laplace transform. Then, with the help of (|8j), we can identify the 



right-hand side of (36) as a sum of $3 functions, which allows us to write (after some algebra) 



4> 3 (6, c; xt, yt) 



l)xt 



1, c - 2; xt, yt) - (c - 2)<l 3 (b - 1, c - 1; xt, yt) - yt<S 3 (6 - 1, c; xt, yt)) (37) 
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which shows that $3 can be recursively expressed via lower values of its first argument. Thus, ( 37 1 can 
be recursively applied to yield 

2(6-1) 

$3 (b,c;xt,yt) = ^ ^$3 (1, c - i; xt, yt) (38) 

i=0 

where on are certain coefficients associated to 63 (1, c— i;xt,yt) that can be obtained by working out 



the recursion (37 1. An explicit formula for 014 can be inferred from the first coefficients on, i = 0, 1, 2, 
leading to 

Finally, ([T0])-([rT]) follow from (f38])-([39l) after setting xt = w, yt = z, and defining the polynomial 



Ai(b,c;z) = (w b ' 1 /z b - 1 - i )a i . 
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